Preparation and spectroscopy of a metastable Mott insulator state with attractive interactions 
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We prepare and study a metastable attractive Mott insulator state formed with bosonic atoms in a three- 
dimensional optical lattice. Starting from a Mott insulator with Cs atoms at weak repulsive interactions, we 
use a magnetic Feshbach resonance to tune the interactions to large attractive values and produce a metastable 
state pinned by attractive interactions with a lifetime on the order of 10 seconds. We probe the (de-)excitation 
spectrum via lattice modulation spectroscopy, measuring the interaction dependence of two- and three-body 
bound state energies. As a result of increased on-site three-body loss we observe resonance broadening and 
suppression of tunneling processes that produce three-body occupation. 



PACS numbers: 37.10.Jk, 67.85.Hj, 03.75.Lm, 05.30.Rt 

Ultracold atomic gases in optical lattices provide a plat- 
form for investigating novel many-body dynamics in a highly 
controllable environment |1|. Recent studies of the quan- 
tum phase transition between a superfiuid and an insulating 
Mott state for bosons [2-6 1 and the metal to insulator transi- 
tion for fermions Q [8) exemplify the control available over 
these systems by tailoring the lattice potential, or tuning inter- 
particle interactions using magnetic Feshbach resonances (9). 
This opens the door towards quantitative studies of phenom- 
ena that are not well understood in condensed matter physics, 
and also novel dynamics beyond what is normally realizable 
in solid state systems. Key examples of the latter include non- 
equilibrium dynamics such as transport processes iflOl fTTI . 
quenches across phase transitions lfl2Ul4 |. and the potential to 
realize metastable many-body states that are long-lived on ex- 
perimental timescales. For example, repulsively bound atom 
pairs can be formed in experiments [TT31 [161 and are stable in 
optical lattice systems because of the lack of strong dissipative 
mechanisms such as lattice phonons that could remove energy 
on short timescales. 

Here we prepare and study a novel metastable many-body 
state, specifically a metastable Mott insulator, in which par- 
ticles are exponentially localized at individual sites through 
attractive two-body interactions. This state can be prepared 
via a sudden change in interparticle interactions, starting from 
a Mott insulator state with repulsive interactions, and then 
switching to attractive interactions abruptly on the tunneling 
timescale in the lattice. This is made possible by a broad Fesh- 
bach resonance for Cs atoms in the lowest hyperfine state ifTTl . 
We demonstrate that the resulting highly excited many-body 
state is long-lived, allowing for detailed studies of its prop- 
erties. Using amplitude modulation (AM) of the optical lat- 
tice corresponding to two-photon Bragg transitions (TH), we 
measure the (de-)excitation spectrum. By identifying specific 
excitation resonances, we map out two-body bound state ener- 
gies over a wide range of scattering lengths and make a quan- 
titative comparison with the corresponding theoretical predic- 
tion Q3D - We also observe shifted resonances connected to 



three-body bound states [20|, which feature a fast broaden- 
ing and strength reduction for increasing attractive interac- 
tions due to increasing three-body loss ETI . This reduction in 
strength corresponds to an inhibition of tunneling events that 
create three-body occupation due to on-site loss processes, 
which promises interesting effects on the many-body physics 
of the system ll22l . 

The standard Bose-Hubbard (BH) model describes the dy- 
namics of bosons in an optical lattice using the two-body on- 
site interaction energy U and the nearest neighbor tunneling 
rate J. There are two groundstates for zero temperature, the 
superfiuid state for U -^iJ and the Mott insulator for U S> J- 
Here we prepare a metastable Mott insulator state with attrac- 
tive interactions (U < 0). This state is analogous to a standard 
bosonic Mott insulator with repulsive interactions, in that par- 
ticles are exponentially localized at different lattice sites, and 
the state exhibits an energy gap due to inter-particle interac- 
tions. However, for U < in a uniform system at unit filling, 
it is the most highly excited Hamiltonian eigenstate of the BH 
model in the limit \U\ J. The state is metastable on long 
timescales because of an energy gap of order U to lower ly- 
ing eigenstates, and there are no fast dissipative mechanisms 
to remove this energy from the system. This is analogous 
to metastability in a gas of repulsively bound pairs |[T31 [T6ll 
and the super- Tonks-Girardeau gas [23|. The preparation of 
this highly excited system is realized by starting in a standard 
Mott insulator with U > 0, U ^ J, and then rapidly switching 
U to a large attractive value, on a timescale that is faster than 
h I J, but sufficiently slow to avoid population of higher Bloch 
bands. 

We start with the production of an essentially pure Bose- 
Einstein condensate (BEC) of 1.0 x 10 5 Cs atoms in an 
optical dipole trap in the lowest hyperfine state following 
mainly the procedures detailed in Ref. 1241 . We load the BEC 
over the course of 500 ms into a cubic optical lattice cre- 
ated by three retro-reflected laser beams with a wavelength of 
A = 1064.5 nm. The lattice depth Vq for the measurements that 
follow is typically set to Vo = 20£'r, where E^ = h 2 /(2mX 2 ) 
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FIG. 1: (color online) Stability of the many-body system for at- 
tractive interactions a) BEC fraction xggc as a function of as with 
/h = 50ms. The narrow shaded areas indicate the locations of Fesh- 
bach resonances. The inset shows the region around as=0afJ. b) 
Number of remaining atoms Na as a function of as with the same 
settings as in a), c) xbec as a function of ?h f° r a s = + 220 ao (cir- 
cles), as = — 240 ao (squares) and as = — 2000 ao (diamonds). The 
solid lines are exponential fits giving 1 /e-decay times of 14.2 ± 0.6 s 
(+220ao), 13.7 ± 0.6s (-240a ) and 11.4±0.4s (-2000a ). d) 
Na as a function of tn for the same settings as in c). The solid 
lines are exponential fits giving decay times of 28. 3 ± 0.9 s (+220 ao), 
25. 2 ± 1.8s (-240 a ) and 19.5 ± 0.5 s (-2000a ). The vertical error 
bars reflect the one-sigma statistical error. 



is the atomic recoil energy with the mass m of the Cs atom. 
The scattering length as can be tuned in a range between 
rj — 2500 ao and ps 1500 ao by applying a magnetic field B be- 
tween 0G and 48 G, employing a broad Feshbach resonance 
with pole at ~ — 12 G and a narrow resonance at ~ 48 G 11251 . 
Here, ao is Bohr's radius. The zero crossing for as is at 
17. 1 19G. During loading of the lattice the interaction strength 
is set to as = +220ao to prepare an atomic Mott insulator state 
as the starting point for all subsequent measurements. The 
number of doubly occupied lattice sites can be controlled via 
the external confinement, which is primarily set by the dipole 
trap laser beams used for the initial BEC preparation. 

To first investigate the stability of the many-body system 
with attractive interactions, we ramp B and therefore the in- 
teraction strength to the desired value with a ramp speed of 
2.5G/ms, wait for a variable hold time fn and return to the 
initial interaction strength. Subsequently we ramp down the 
optical lattice and recapture the cloud of particles in the dipole 
trap. Using the time-of-flight technique, we deduce the num- 
ber of remaining atoms Na and the BEC fraction xbec- ln the 
course of the ramps, we have to cross several narrow Feshbach 
resonances [17|. To avoid heating through interaction-induced 



FIG. 2: (color online) Excitation spectrum for a Mott insulating 
state with attractive interactions a) BEC fraction xg^c as a func- 
tion of /m for high (low) initial density shown as circles (squares) 
at a$= — 306 ao- The solid lines are double Gaussian fits, b) Nor- 
malized molecule number Nm as a function of /m at as = — 306 ao. 
The vertical error bars reflect the one-sigma statistical error, c) Set 
of excitation spectra with x^ec as a function of /m as as is varied: 
as = — 1 1 1 ao (top) to as = — 2334 ao (bottom). The three-body res- 
onance is indicated by an arrow. 



band-transfer [26 1, we cross the resonances using fast ramps 
with a speed of 2 x 10 4 G/ms l27l . Figure[IJa) and (b) show 
xbec and Na as a function of as for fH=50ms. Sharp dips 
for both observables occur near the expected locations of the 
narrow resonances. Also, around the zero crossing of as the 
system becomes unstable. An approximately 20% (10%) de- 
crease for xbec C^a) is observed for as <0 compared to the 
values at as >0, probably caused by the crossing of the zero 
interaction region within a finite time. Figure[TJc) and (d) 
show xbec and Na as a function of fn for different values of 
as- As one would expect, xbec decays somewhat faster than 
Na at given as. Exponential fits to the data yield lifetimes. 
These are reduced for more attractive interactions, but stay in 
the range of 10 to 20 seconds, more than 30 times larger than 
the bare tunneling time h/Jai this lattice depth. This clearly 
shows that the system is stabilized by attractive interactions 
and allows us to investigate its properties in more detail. 

In this work we focus on the (de-)excitation spectrum 
12 [T8l . which we measure by AM at frequency /m of one 
of the lattice beams at typically 20 % of Vq for a duration of 
fH = 300ms. Tunneling processes at U < to sites with non- 
zero occupation lower the overall energy and lead to a deex- 
citation of the system. The deexcitations - through the spatial 
configuration - are mapped onto excitations when returning 
subsequently to repulsive interactions, leading to an increased 
overall energy of the system, which we detect again by mea- 
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suring xbec as a sensitive indicator for the energy deposited 
into the system. Figure[2|a) shows the measured excitation 
spectrum for as = — 306 ao in the vicinity of \U\/h « 1.6kHz 
as calculated from the standard BH model. A double reso- 
nance structure can clearly be identified, similar to our previ- 
ous work in the regime of strong repulsive interactions 11281 . 
At comparatively low initial densities, with a small number 
of doubly occupied sites, the lower resonance is dominant, 
whereas for higher initial densities, giving a larger fraction 
of doubly occupied sites, the upper resonance becomes more 
pronounced at the expense of the lower resonance. As in the 
repulsive case, this splitting is caused by the energy differ- 
ence for excitations in the different Mott shells through ef- 
fective multibody interactions for three particles at the same 
lattice site [20]. We associate the lower resonance with exci- 
tations in the singly occupied shell, creating doubly occupied 
sites, i. e. two-body bound states in the presence of the lat- 
tice, and the upper resonance with excitations in the doubly 
occupied shell, creating triply occupied sites, i. e. three-body 
bound states in the presence of the lattice. In contrast to the 
repulsive case {28 1 the three-body resonance is now at higher 
energies. Our interpretation is confirmed by a measurement 
of the number of doubly occupied sites through molecule for- 
mation l27l . Figure|2|b) shows the molecule number A^m as 
a function of /m. For the lower resonance we find a greatly 
increased probability for dimer formation, which we attribute 
to an increase in the number of sites with double occupancy, 
while for the upper resonance we observe a strongly reduced 
probability for dimer formation, in agreement with the fact 
that doubly occupied sites are resonantly emptied upon mod- 
ulation and that particles at triply occupied sites are lost due 
to fast three-body recombination at negative as ET1 . 

Figure|2|c) shows excitation spectra taken over a compara- 
tively large frequency range as as is varied from — 112ao to 
— 2334 ao. The resonance discussed above corresponding to 
the excitation into the two-body bound state is clearly visible 
and it is shifted to higher frequencies as |«s| is increased, as 
one would expect. The three-body resonance is only visible 
for weak interaction strengths. It also shifts to higher fre- 
quencies, as expected. We will discuss the behavior of this 
resonance below. Interestingly, we observe another resonance 
with an inverse behavior compared to the two-body and three- 
body resonances, with decreasing frequency as \as \ increases. 
We observe an additional resonance at twice the frequency 
with the same behavior, visible in the spectra taken for strong 
interactions. These resonances do not disappear when we pre- 
pare a purely singly-occupied Mott insulator. We therefore 
suspect that they are related to excitations into the first excited 
two-body bound state of the lattice [19|. The strong visibility 
of the half-frequency resonance most likely is a result of the 
presence of the first harmonic in the frequency spectrum of J, 
which for our parameters reaches a strength of up to 20 % of 
the main frequency component 11291 . 

Figure[3] provides an overview over the measured excita- 
tion energies AE as as is varied and compares our data to the 
prediction by theory |[T9l 1281 . The three relevant excitation 
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FIG. 3: (color online) Excitation energies AE into two- and three- 
body bound states. Two-body bound state excitation energy (circles), 
three-body bound state excitation energy (squares), and first excited 
two-body bound state excitation energy (diamonds, obtained by dou- 
bling the measured values) are shown as a function of ag. We have 
added to the plot the data for ag > from Ref. 1281 . augmented by 
new measurements for as > 1000aQ. The shaded areas at the top 
indicate the first two single-particle bands. The solid (dashed) line 
is the calculated excitation energies for the two-body (three-body) 
bound state and the dotted line is the calculated excitation energy for 
the first excited two-body bound state. The dashed-dotted line is the 
standard BH model calculation. The vertical error bars reflect the 
one-sigma statistical error as derived from the gaussian fit. The hori- 
zontal error bars indicate the variation of as over the cloud due to the 
gradient in B. The inset illustrates the main excitation processes. 



processes are illustrated in the inset to this figure. We deter- 
mine AE — Ii/m.c by fitting simple gaussians to the loss fea- 
tures as shown in Figure[2ja) and taking the peak positions 
/m,c- The data for repulsive interactions (as > 0) is taken 
from Ref. (28), augmented by new measurements at strong re- 
pulsive interactions (as > 1000 ao). In general, we find good 
agreement between our measurements and the calculated en- 
ergies for the excitations to the lowest-band two-body bound 
state using the exact numerical results [ 19 29 1, and the fitting 
function described in Ref. [28], which extrapolates the result 
for the three-body resonance from Ref. Ifl9l . Only for compar- 
atively strongly attractive interactions do we find a significant 
deviation for the energy of the three-body bound state, which 
is expected due to more complex three-body physics arising 
at negative values for the scattering length [ 30 1 . The exci- 
tation energies into the first excited two-body bound states, 
derived by doubling the measured frequency values, show the 
same qualitative behavior as the exact numerical calculations 
|fT9"l|29 |, though with a significant offset for strong attractive 
interactions. Note that the exact numerical analysis is only 
valid for two particles or a very dilute system, and does not 
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FIG. 4: (color online) Broadening of the three-body excitation res- 
onance, a) Set of excitation spectra showing the normalized BEC 
fraction xbec as a function of /m for various values of as as indi- 
cated. The resonance peak connected to the two-body (three-body) 
bound state is located at lower (higher) frequency. The solid lines 
are double-peaked Gaussian fits, from which the width (FWHM) of 
the resonances is deduced, b) FWHM of the three-body excitation 
resonance as a function of as for attractive interactions, c) For com- 
parison, we show the FWHM of the three-body excitation resonance 
as a function of as for repulsive interactions [28|. 



include the influence of a Mott insulating background. Its 
contribution will increase for stronger particle interactions and 
might account for the offset observed. Also, we are not able 
to completely exclude the possibility of a systematic deviation 
in the calculation of as in this magnetic field region. 

Detailed spectra with the three-body bound state excitation 
resonance for attractive interactions are shown in Fig.|4ja). 
Compared to the two-body bound-state resonance at lower fre- 
quencies, we clearly observe a broadening of the three-body 
resonance as as is increased. In fact, the resonance broadens 
so much that it finally disappears. Given our finite signal- 
to-noise ratio, we are not able to detect the resonance for 
|a s | > 500a . Fig.^b) shows the width (FWHM) as a func- 
tion of as. We note that, as there is currently no theory on the 
lineshapes of three-body resonances, the FWHM is obtained 
from simple gaussian fits as used above. The broadening leads 
to a decrease in the height of the resonance, i.e., a decrease in 
the formation rate of triply-occupied sites, and thus to a de- 
crease in the overall loss rate. We interpret this behavior in 
terms of the process discussed in Ref. Il22ll : Fast three-body 
loss can suppress the formation of triply occupied sites, anal- 
ogous to the suppression of inelastic two-body processes as 
reported in Ref. [31] for the case of weakly-bound molecules. 
As the three-body recombination rate scales like N ~ a^n 3 
with the density «, one would expect a strong increase of the 



three-body loss rate. We note that for repulsive interactions 
a similar broadening of the three-body bound-state excitation 
resonance can be observed. This is shown in Figure|4|c). In 
this case the increase of the FWHM as a function of as is far 
less drastic. Further experimental and theoretical investiga- 
tions will be needed for a quantitative analysis of the broaden- 
ing, including a calculation of the modified on-site density of 
three-body bound states and the corresponding tunneling rates 
at attractive and repulsive interactions and the role of Efimov 
physics in confined dimensions [32 1. 

In summary we have investigated a metastable Mott insulat- 
ing state with attractive interactions and have shown that this 
state exhibits a very long lifetime on the order of 10 seconds. 
By measuring excitation spectra we were able to determine 
the energies of two-body and three-body bound states. The 
broadening of the excitation resonance corresponding to the 
three-body bound state gives a strong indication that very high 
three-body recombination rates suppress the creation of triply 
occupied sites, thereby reducing the effective loss rate of the 
system, which also can be interpreted in terms of the quan- 
tum Zeno effect [ 33 1 with the three-body recombination rate 
playing the role of the measurement. This effect can be used 
to realize a Bose-Hubbard model with a three-body hard-core 
constraint [22], leading to interesting many-body physics, in- 
cluding a dimer superfluid phase and a continuous supersolid 
phase for attractive bosons l34ll35ll . the realization of Pfaffian- 
like states in one-dimensional geometries [36|, and the stabi- 
lization of an atomic color superfluid for fermions ||37l . The 
Mott insulating state with attractive interactions can also serve 
as starting point for the preparation of stable superfluid con- 
densates at finite momentum with negative temperature [38]. 

We are indebted to R. Grimm for generous support. We 
thank J. von Stecher, P. Johnson, and E. Tiesinga for fruitful 
discussions. We gratefully acknowledge funding by the Aus- 
trian Science Fund (FWF) within project I153-N16 and within 
the framework of the European Science Foundation (ESF) Eu- 
roQUASAR collective research project QuDeGPM. 



[1] I. Bloch, J. Dalibard, and W. Zwerger, Rev. Mod. Phys. 80, 885 
(2008). 

[2] M. Greiner, O. Mandel, T. Esslinger, T. W. Hansch, and I. Bloch, 

Nature 415, 39 (2002). 
[3] N. Gemelke, X. Zhang, C.-L. Hung, and C. Chin, Nature 460, 

995 (2009). 

[4] T. Fukuhara, S. Sugawa, M. Sugimoto, S. Taie, and Y. Takahashi, 

Phys. Rev. A 79, 041604 (2009). 
[5] W. S. Bakr, A. Peng, M. E. Tai, R. Ma, J. Simon, J. I. Gillen, S. 

Foiling, L. Pollet, and M. Greiner, Science 329, 547 (2010). 
[6] J. F. Sherson, C. Weitenberg, M. Endres, M. Cheneau, I. Bloch, 

and S. Kuhr, Nature 467, 68 (2010). 
[7] R. Jordens, N. Strohmaier, K. Giinter, H. Moritz, and T. 

Esslinger, Nature 455, 204 (2008). 
[8] U. Schneider, L. Hackermuller, S. Will, Th. Best, I. Bloch, T. A. 

Costi, R. W. Helmes, D. Rasch, and A. Rosch, Science 322, 1520 

(2008). 



[9] C. Chin, R. Grimm, P. S. Julienne, and E. Tiesinga, Rev. Mod. 

Phys. 82, 1225 (2010). 
[10] L. Hackermiiller, U. Schneider, M. Moreno-Cardoner, T. Kita- 

gawa, T. Best, S. Will, E. Dernier, E. Altman, I. Bloch, and B. 

Paredes, Science 327, 1621 (2010). 
[11] S. Trotzky, Y.-A. Chen, A. Flesch, I. P. McCulloch, U. 

Schollwock, J. Eisert, and I. Bloch (201 l), |arXiv:1101.2659| 
[12] M. Cheneau, P. Barmettler, D. Poletti, M. Endres, P. SchauB, 

T. Fukuhara, C. Gross, I. Bloch, C. Kollath, and S. Kuhr (2011), 

larXiv:1111.0776H . 
[13] D. Chen. M. White. C. Borries. and B. DcMarco. Phys. Rev. 

Lett. 106, 235304(2011). 
[14] C. Kollath, A. M. Liiuchli, and E. Altman, Phys. Rev. Lett. 98, 

180601 (2007). 

[15] K. Winkler, G. Thalhammer, F. Lang, R. Grimm, J. Hecker 

Denschlag, A. J. Daley, A. Kantian, H. P. Buchler, and P. Zoller, 

Nature 441, 853 (2006). 
[16] N. Strohmaier, D. Greif, R. Jordens, L. Tarruell, H. Moritz, T. 

Esslinger, R. Sensarma, D. Pekker, E. Altman, and E. Demler, 

Phys. Rev Lett. 104, 080401 (2010). 
[17] C. Chin, V. Vuletic, A. J. Kerman, S. Chu, E. Tiesinga, P. J. Leo, 

and C. J. Williams, Phys. Rev. A 70, 032701 (2004). 
[18] T. Stoferle, H. Moritz, C. Schori, M. Kohl, and T. Esslinger, 

Phys. Rev. Lett. 92, 130403 (2004). 
[19] H. P. Buchler, Phys. Rev Lett. 104, 090402 (2010). 
[20] P. R. Johnson, E. Tiesinga, J. V. Porto, and C. J. Williams, New 

J. Phys. 11, 093022 (2009). 
[21] T. Kraemer, M. Mark, P. Waldburger, J. G. Danzl, C. Chin, B. 

Engeser, A. D. Lange, K. Pilch, A. Jaakkola, H.-C. Nagerl, and 

R. Grimm, Nature 440, 3 1 5 (2006) . 
[22] A. J. Daley, J. M. Taylor, S. Diehl, M. Baranov, and P. Zoller, 

Phys. Rev. Lett. 102, 040402 (2009). 
[23] E. Haller, M. Gustavsson, M. J. Mark, J. G. Danzl, R. Hart, G. 

Pupillo, and H.-C. Nagerl, Science 325, 1224-1227 (2009). 
[24] T. Kraemer, J. Herbig, M. Mark, T. Weber, C. Chin, H.-C. 

Nagerl, and R. Grimm, Appl. Phys. B 79, 1013 (2004). 
[25] A. D. Lange, K. Pilch, A. Prantner, F. Ferlaino, B. Engeser, 

H.-C. Nagerl, R. Grimm, and C. Chin, Phys. Rev. A 79, 013622 

(2009). 

[26] M. Kohl, H. Moritz, T. Stoferle, K. Gunter, and T. Esslinger, 

Phys. Rev. Lett. 94, 080403 (2005). 
[27] J. G. Danzl, M. J. Mark, E. Haller, M. Gustavsson, R. Hart, A. 

Liem, H. Zellmer, and H.-C. Nagerl, New J. Phys. 11, 055036 

(2009) . 

[28] M. J. Mark, E. Haller, K. Lauber, J. G. Danzl, A. J. Daley, and 
H.-C. Nagerl, Phys. Rev. Lett. 107, 175301 (201 1). 

[29] See Supplemental Material at [URL will be inserted by pub- 
lisher] 

[30] J. von Stecher, private communication (2011). 

[31] N. Syassen, D. M. Bauer, M. Lettner, T. Volz, D. Dietze, J. J. 

Garc ia-Ripoll, J. I. Cirac, G. Rempe, and S. Diirr, Science 320, 

1329 (2008). 

[32] J. Portegies and S. Kokkelmans, Few-Body Systems 51, 219 
(2011). 

[33] B. Misra and E. C. G. Sudarshan, J. Math. Phys. 18, 756 (1977). 
[34] S. Diehl, M. Baranov, A. J. Daley, and P. Zoller, Phys. Rev. Lett. 

104, 165301 (2010). 
[35] L. Bonnes and S. Wessel, Phys. Rev. Lett. 106, 185302 (2011). 
[36] B. Paredes, T. Keilmann, and J. I. Cirac, Phys. Rev. A 75, 

053611 (2007). 

[37] A. Kantian, M. Dalmonte, S. Diehl, W. Hofstetter, P. Zoller, and 

A. J. Daley, Phys. Rev. Lett. 103, 240401 (2009). 
[38] A. Rapp, S. Mandt, and A. Rosch, Phys. Rev. Lett. 105, 220405 

(2010) . 



6 



SUPPLEMENTARY MATERIAL: PREPARATION AND 
SPECTROSCOPY OF A METASTABLE MOTT INSULATOR 
STATE WITH ATTRACTIVE INTERACTIONS 

Modifications and frequency components of U and J during 
amplitude modulation 

A sinusoidal modulation of the lattice depth leads to a mod- 
ulation in both U and J. Excitations are driven mainly by 
the modulation of J 0"). To characterize this process in a 
first approximation, we write the lattice depth modulation as 
V(t) = Vb(l +Msin(o)Mf )) with 0)m = 2tt/m and numerically 
calculate U (?) and J{t) for each time step using tunnelling 
rates and Wannier functions for the lowest Bloch band. Here, 
we do not include the modification of U and J arising for 
stronger interactions. To derive the frequency components of 
7(f), we Fourier transform J(t), omit the DC part and nor- 
malize the strength of the frequency components to the peak 
strength at the main frequency o>m- Due to the nonlinear 
dependence of J on Vb, higher harmonics at 2©m and 3o)m 
are visible, and their strength as a function of the modula- 
tion depth M is shown in figure[5|a) for three different val- 
ues of the lattice depth Vq. For our experimental parameters 
(Vb = 20 E R , M = 20%) the 2c^m component reaches 20% 
peak strength of the main frequency peak, whereas the 3o>m 
component stays below 3%. This explains the good visibil- 
ity of the half-frequency resonance into the first excited two- 
body bound state. Note that a half-frequency resonance is 
also present for the lowest two-body bound state as shown in 
Ref. 0. The weak 3©m component would give rise to a fur- 
ther resonance, which, however, is not observable given the 
finite signal-to-noise ratio for our experiment. 




Modulation depth M (%) Modulation depth M (%) 

FIG. 5: (color online) a) Normalized strength of the 2(0 and 3co fre- 
quency components of J(t) as a function of the modulation depth M 
for Vb = 20Er (solid line), Vo = 15Er (dashed line), and Vb = 25 Er 
(dotted line), b) Normalized change for the time averaged values (U) 
and (J) for U and J as a function of M for different values of Vo as in 
a). 

Due to the nonlinearity of J and U as a function of Vb the 
corresponding time averaged values during the modulation are 
different from the values that one obtains assuming no modu- 
lation. Figure[5|b) shows the relative change of the time aver- 
aged (U) and (J) as a function of M for three different values 
of Vq. The shift of U by less than 0.5 % is nearly independent 
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cut-off parameter A 

FIG. 6: (color online) Convergence of the numerical analysis for in- 
creasing cut-off parameter and number of Bloch bands: Regularized 
shift of the detuning Xh{E) at energy E = 5.5£r and for an optical 
lattice with Vo = 16£r. The dots are the values for different shell pa- 
rameters 5. The shift satisfies the asymptotic behavior Xx = C /A + X 
(solid line). The inset shows the fast convergence for fixed A by in- 
creasing the shell parameter S, i.e., the involved number of Bloch 
bands. 

of Vo and can be neglected in view of the uncertainties in our 
experiment. The shift of J is more pronounced, but it does not 
affect the experimental results on the measurement of inter- 
action energies. Nonetheless, this effect has to be taken into 
account for a future investigation of modulation-assisted tun- 
neling rates, additionally to the changes due to modifications 
in the on-site wave functions arising from interactions [3]. 



Numerical analysis of the two-particle bound states 

The numerically efficient and precise determination of the 
attractive and repulsive bound state energies for two particles 
in a three-dimensional optical lattice has been presented in 
Ref. J?) . Here, we use this method for the comparison with 
the experimental data; the notation is in analogy to Ref. Pfl . 
Within the two channel description of the short range pseudo- 
potential, the s-wave scattering length «s 1S related to the 
detuning v of the molecular channel, and the coupling g 
via 47th 2 as/m = —g 2 /v. The influence of the optical lat- 
tice is a shift %(E) in the detuning and is associated with a 
change of the free particle properties due to the formation 
of Bloch bands. The relation between the bound state en- 
ergies E bs and the shift %(E) in the detuning takes the form 
%(E bk ) — —7ta/(8Enas) with a the lattice spacing. The nu- 
merical determination of %(E) involves the precise determi- 
nation of Bloch wave functions and band energies, as well 
as a summation over a high number of Bloch bands. In ad- 
dition, it requires a regularization of the coupling between 
the open and closed channel. Here, we choose the regu- 
larization a(r) = J v ( A jc/kexp(/kr)/(27i:) 3 , where the volume 

v(A) = A 3 vo is centered around k = with vo = (2ft) 3 /a 3 
the volume of the first Brillouin zone. For fixed cut-off A, 
the summation over Bloch bands converges very quickly for 
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S > A, see inset to Fig. |6j the shell parameter S denotes the 
number of Bloch bands for each spatial direction included 
in the summation. Finally, we can remove the cut-off A via 
the asymptotic scaling relation Xa{E) = c / A + %{E) (5), see 
Fig. [6] The convergence of the numerical data is checked by 
varying the shell parameter S, cut-off A, and the precision 
for the determination of the Bloch wave functions and Bloch 
bands. 

Using the above procedure, the repulsive and attractive 
bound states are determined for the experimentally relevant 
parameters with Vb = 20£r. The lowest bound state for weak 
interactions is dominated by the contributions of two parti- 
cles in the lowest Bloch band, and for increasing interactions 
smoothly connects to the state with a single bound molecule 
of mass 2m within the lowest Bloch band of an optical lat- 
tice with strength 2Vo. In turn, the first excited bound state 
corresponds for weak interactions to a state with one parti- 



cle in the lowest Bloch band and the second particle in the 
first excited Bloch band. Again for increasing interactions, 
this state smoothly connects to a single molecule in the first 
excited Bloch band. 
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